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1 Introduction 

Electron-positron Bhabha scattering plays a special role in particle phenomenology. It pro- 
vides a very efficient tool for luminosity determination at electron-positron colliders and 
thus it is crucial for extracting physics from the experimental data. Small angle Bhabha 
scattering has been particularly effective as a luminosity monitor at the energies of LEP and 
SLC because its cross section is large and QED dominated |l|2j . At a future International 
Linear Collider (ILC) the luminosity spectrum is not monochromatic due to beam-beam 
effects. Therefore measuring the cross section of the small angle Bhabha scattering alone 
is not sufficient, and the acollinearity of the large angle Bhabha scattering has been sug- 
gested for disentangling the luminosity spectrum |3|4j . Large angle Bhabha scattering is 
important also at colliders operating at a center of mass energy of a few GeV, such as 
BABAR/PEP-II, BELLE/KEKB, BES/BEPC, KLOE/DA$NE, and VEPP-2M, where it is 
used to measure the luminosity j^]. Since the accuracy of the theoretical evaluation of the 
Bhabha cross section directly affects the luminosity determination, remarkable efforts have 
been devoted to the study of the radiative corrections to this process (see [T] for an extensive 



1 



list of references). Pure QED contributions are particularly important because they domi- 
nate the radiative corrections to the large angle scattering at intermediate energies 1-10 GeV 
and to the small angle scattering also at higher energies. The calculation of the QED radia- 
tive corrections to the Bhabha cross section is among the classical problems of perturbative 
quantum field theory with a long history. The first order corrections are well known (see 
|6|7j and references therein). To match the impressive experimental accuracy the complete 
second order QED effects have to be included on the theoretical side. The evaluation of 
the two-loop virtual corrections constitutes the main problem of the second order analysis. 
The complete two-loop virtual corrections to the scattering amplitudes in the massless elec- 
tron approximation have been computed in Ref. [B], where dimensional regularization has 
been used for the infrared divergences. However, this approximation is not sufficient since 
one has to keep a nonvanishing electron mass to make the result compatible with available 
Monte Carlo event generators |l|5|9|lU|llj . Recently an important class of the second order 
corrections, which include one closed fermion loop, has been obtained for a finite electron 
mass ^2] including the soft photon bremsstrahlung ^3]- A similar evaluation of the purely 
photonic two-loop corrections is a challenging problem at the limit of present computational 
techniques |14|15|16j . The most complete result available so far can be found in Ref. [T7] 
where the contribution of double box diagrams is still missing. On the other hand in the 
energy range under consideration only the leading contribution in the small ratio m? e / s is of 
phenomenological relevance and should be retained in the theoretical evaluations. For arbi- 
trary scattering angle even in this approximation only the two-loop corrections enhanced by 
a power of the large logarithm ln(mg/s) are known so far |18|19j . In the limit of the small 
scattering angle, however, the structure of the corrections is much simpler [20 that allowed 
for the evaluation of the corrections up to the nonlogarithmic term |21|22j . The result for 
the nonlogarithmic contribution for arbitrary scattering angle has been reported in a letter 
23 . It was obtained by employing the general theory of infrared singularities in QED which 
allows to reduce the calculation in the small electron mass approximation to the analysis of 
a strictly massless scattering amplitude and the massive vector form factor. In the present 
paper we describe the details of this calculation. In the next section we outline the structure 
of the perturbative expansion for the Bhabha cross section. In Sect.|3]we consider the struc- 
ture of the infrared logarithms and formulate the method of infrared subtractions. In Sect. 0] 
the explicit relation between the amplitudes of the massive and massless Bhabha scatter- 
ing is established through the infrared matching procedure and the result for the two-loop 
corrections to the massive Bhabha scattering is obtained. Sect. |S] contains the numerical 
estimates and the summary. 

2 Perturbative expansion of the cross section 

We consider the phenomenologically interesting kinematical region s, t, m> m 2 el where all 
the terms suppressed by the electron mass can be neglected. The perturbative expansion for 



2 



the Bhabha cross section in the fine structure constant a is defined as follows 
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In the small electron mass approximation the leading order differential cross section takes 
the form 

do-W a 2 /l-x + x 2 ^ 2 
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where x = (1 — cos6 l )/2 and # is the scattering angle. 

The virtual corrections taken separately suffer from the so/t divergences, which can be 
regulated e.g. by giving the photon a small auxiliary mass A. These soft divergences are 
canceled in the inclusive cross section when one adds the photonic bremsstrahlung |24|25j . 
The standard approach to deal with the bremsstrahlung is to split it into a soft part which 
accounts for the emission of the photons with the energy below some cutoff e cut <C m e , 
and a hard part corresponding to the emission of the photons with the energy above e cu t. 
The infrared finite hard part is then computed numerically using Monte-Carlo methods 
with physical cuts dictated by the experimental setup. At the same time the soft part is 
computed analytically and combined with the virtual corrections ensuring the cancellation of 
the singular dependence on A. Note that in many practical realizations of the Monte-Carlo 
event generators the cancellation of the infrared singularities is build in and implemented to 
high orders of perturbation theory for the amplitudes rather than for the cross section (see 
e.g. |9|lUj ). We will come back to this issue in Sect. 14.21 Thus in the first order we consider 
the sum of one-loop virtual correction and single soft photon emission 

= #> + 50) , (3) 

where 
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and the expressions for 80) and 80) are given by Eqs. 1)4*51 fl""")) of the Appendix. Eq. (JSJ) can 

be decomposed according to the asymptotic dependence on the electron mass 
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For the pure photonic correction the coefficients 8^ read (see e.g. |18pi9j ) 
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Li n (z) is the polylogarithm, e = y/s/2, e cut is the energy cut on the emitted soft photon. 
The second order correction can be represented as a sum of three terms 

S {2) = <© + S& + S$ (8) 

which correspond to the two-loop virtual correction including the one-loop corrections to the 
amplitude square, one-loop virtual correction to single soft photon emission, and the double 
soft photon emission, respectively. In the small electron mass limit it has the following 
decomposition 
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The pure photonic, i. e. without closed fermion loops, logarithmically enhanced contribution 
reads 
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where C(3) = 1.202057. . . is the value of the Riemann's zeta-function. In the rest of the 

(2) 

paper we focus on the photonic contribution to 8q . 



3 Structure of infrared logarithms 

The general problem of the calculation of the small electron mass asymptotics of the correc- 
tions including the power-suppressed terms can systematically be solved within the expansion 
by regions approach [26 ;27j- We, however, are interested only in the leading order term. The 
leading order contribution in Eq. (JSJ) contains the logarithmic terms, which become singular 
as m e approaches zero revealing the collinear divergences regulated by the electron mass. In 
the massless limit both the collinear and the soft divergences can be treated by dimensional 
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regularization as well. Here we should note that the collinear divergences in the massless 
approximation are also canceled in a cross section which is inclusive with respect to real 
photons and electron-positron pairs collinear to the initial or final state fermions |28j . This 
means that if an angular cut on the collinear emission is sufficiently large, 6 cut ^> \Jm 2 e / 's, the 
inclusive cross section is insensitive to the electron mass and can in principle be computed 
with m e = by using dimensional regularization for the infrared divergences for both virtual 
and real radiative corrections like it is done in the theory of QCD jets. However, as it has 
been mentioned above, all the available Monte Carlo event generators for Bhabha scattering 
with specific cuts on the photon bremsstrahlung dictated by the experimental setup employ 
a nonzero electron mass as an infrared regulator, which therefore has to be used also in the 
calculation of the virtual corrections. As far as the leading term in the small electron mass 
expansion is considered, the difference between the massive and the dimensionally regular- 
ized massless Bhabha scattering can be viewed as a difference between two regularization 
schemes for the infrared divergences. With the known massless two-loop result at hand, the 
calculation of the massive one is reduced to constructing the infrared matching term which 
relates two above regularization schemes. To perform the matching we develop the method 
of infrared subtractions which simplifies the calculation by fully exploiting the information 
on the general structure of infrared singularities in QED. The method was originally ap- 
plied in Ref. to the analysis of the two-loop corrections to the vector form factor in 
an Abelian gauge model with mass gap. Let A^ 2 \m ei X) be the two-loop contribution to 
the massive electron-positron scattering amplitude with the photon mass used to regulate 
the soft divergences. The main idea of the method is to construct an auxiliary amplitude 
*4( 2 )(m e , A), which has the same structure of the infrared singularities but is sufficiently sim- 
ple to be evaluated at least in leading order in the small mass expansion. Then the difference 
4(2) _ ^(2) hag a finite limit £4(2) 

as m e , A tend to zero. This quantity does not depend on 
the regularization scheme for A^ and A^ 2 \ It and can be evaluated by using dimensional 
regularization for each term and then taking the limit of four space-time dimensions. The 
full amplitude is given by a sum 

A i2) (m e , A) = A {2 \m e , A) + 5A {2) + 0(m e , A) . (11) 

Thus the infrared divergences, which induce the asymptotic dependence of the virtual cor- 
rections on the electron and photon masses, are absorbed into the auxiliary amplitude while 
the technically most nontrivial calculation of the term 5A^ is performed in the massless 
approximation. The matching of the massive and massless results is necessary only for the 
singular auxiliary amplitude. Note that the method does not require a diagram-by-diagram 
subtraction of the infrared divergences since only a general information on the infrared struc- 
ture of the total two-loop correction is necessary to construct A^ 2 '(m e ,X). Our analysis is 
based on the following infrared properties of the corrections to the scattering amplitudes: 

(i) exponentiation of the infrared logarithms |3U|31|32|33|34|35| : 

(ii) factorization of the collinear logarithms into external legs j36j ; 
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(iii) nonrenormalization of the infrared exponents |31|33|34] . 



The first two properties are general and hold also for the closed fermion loop corrections 
and non-Abelian gauge theories. The last property is characteristic for the pure photonic 
corrections and plays a crucial role in our analysis. In Sects. I3.1H3.2I by means of (i)-(iii) 
we show that A^(m e , A) for the photonic contributions can be constructed of the two-loop 
corrections to the vector form factor and products of the one-loop corrections. 



3.1 Vector form factor 



The vector form factor T determines the electron scattering amplitude in an external field. It 
plays a special role since it is the simplest quantity which includes the complete information 
about the collinear logarithms, which is directly applicable to a process with an arbitrary 
number of electrons/positrons. Let us consider three different Sudakov asymptotic regimes: 
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A = 0: 



(b) \Q\ > m e > A; 

(c) \Q\ > A > m e ; 

where Q is the Euclidean momentum transfer. In case (a) the soft and collinear divergences 
are treated by dimensional regularization. In case (b) the collinear and soft divergences 
are regularized by m e and A, respectively. In case (c) the photon mass regulates both soft 
and collinear divergences. Though (c) has no direct application to QED, it is instructive 
to study yet another regularization scheme to get deeper insight into the general structure 
of infrared logarithms. We define the perturbative series for the form factor as follows: 
F — J2^=o (f ) ■ The one-loop coefficients read 
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The asymptotic dependence of the form factor on Q in the Sudakov limit is governed by the 
evolution equation [33 34 35 which for the pure photonic contribution takes the form 
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where the anomalous dimension (j)(m e , A, e, a) is a series in a with the coefficients depending 
on the infrared regulators. We can write down the solution of Eq. (JT5|l in the above three 
cases 
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where 0(a) indicates the presence of all order corrections to the coefficients starting with 
O (a) term. In derivation of Eq. (JTSJl we have taken into account that in the case (b) also 
the logarithms of the photon mass exponentiate [31] 
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The exponentiation of the "Sudakov" logarithms is a general property valid also for the 
corrections due to the fermion loops and for non-Abelian gauge theories. The exponent for 
the pure photonic corrections, however, has two distinguished properties. First, the double 
logarithmic term in the exponent is protected against the perturbative corrections. This fact 
is well known since the pioneering works |3U|31l33l34j . A new observation is that beyond the 
first order in a the coefficients of the series for the single logarithmic term in the exponent are 
mass-independent and, therefore, should be the same in all three cases under consideration. 
The derivation of this result to a large extent repeats the proof of the nonrenormalization 
of the double logarithmic contribution. We refrain from giving the details of the derivation 
since it is not directly related to the subject of the present paper. We should emphasize that 
the above properties are not valid for the nonphotonic corrections, e.g. the second order 
single logarithmic contribution of the closed fermion loop to the exponent depends on the 
mass ratio X/m e and is different for the cases (b) and (c) 37J. 

Thus, as far as the photonic corrections are concerned, the perturbative series for the 
logarithm of the form factor beyond one loop includes only the first power of the logarithm 
with the universal coefficients. This means that the coefficients in the series for the form 
factor have the following structure 
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where each C^ n \ n > 1 is equal for (a), (b), and (c). This prediction can be confronted with 
the explicit results for the two-loop corrections which are available in all three cases and read 
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-^ 2 ln 2 (2) + ^ln 4 (2) + 16Li 4 Q . (22) 

In the massless approximation the two-loop correction ()20j) has been known for a long time 
|38l39j . The two-loop correction (j2T|) was first obtained in Ref . |IU] by integrating the disper- 
sion relation with the spectral density computed in Ref. The result has been checked 
in Ref. [23] and can also be found in Ref. [12] as a specific limit of the result for an arbitrary 
momentum transfer. The two- loop correction ()22|) has been obtained in Ref. [2H] • As we see 
the two-loop corrections indeed have the universal logarithmic term corresponding to 

c<2, -!-y + I« 3 >- < 23 > 



By using the recent three-loop result [43J for the massless case we can completely predict 
the three-loop logarithmic corrections in massive cases, which are given by Eq. (fT9*|) with 

C(3) = ~- 3tt 2 - 68C(3) + ^vr 2 C(3) + 240C(5) . (24) 
3.2 Scattering amplitude 

In the high energy limit the amplitude for the electron-positron scattering has two compo- 
nents corresponding to the scattering of particles of the same or opposite chirality. We can 
write the perturbative series for the amplitude as follows 

•*=£(-) ^ (n) ' A^^A^A^^^A^Afl 4 (0) = 1, (25) 

n=0 ^ n ' i=l 

where A^ is a two component vector in the chiral basis corresponding to the tree amplitude. 
The collinear divergences are completely determined by external legs [36] and, therefore, are 
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the same for the scattering amplitude and the square of the form factor. It is convenient to 
introduce a reduced amplitude 

A = T 2 A , (26) 

which is free of collinear divergences. It satisfies a linear differential equation |35|44j which 
for the photonic contribution takes the following form 

9 vU-^lnf-^U (27) 
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where the angular dependent anomalous dimension does not depend on chirality. The solu- 
tion of Eq. (J27)l reads 



•A|a=o = (*4 (0) + 0(a)) exp 
A\wo = (U (0) + 0(a)) exp 




(21 



where the corrections in the prefactor of the exponent are different for different chiral com- 
ponents of the amplitude. There are no photonic corrections to the exponent in Eq. (|28j) 
[3*T] . Note that in the case (b) all the singular dependence of the corrections to the scattering 
amplitude on m e is absorbed into the form factor. Eq. (J28)) implies that the logarithm of 
the reduced amplitude is finite beyond one loop and the coefficients of the series for each 
component of the chiral basis have the following structure 
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Now we can predict the singilar structure of the photonic corrections to the full amplitude 
to all orders and to construct the auxilary amplitude A?>. In the two- loop approximation 
by using Eqs. (Q l2U|) we obtain 
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The expression ()31|) has all the necessary properties: it is composed of the one-loop cor- 
rections to the chiral amplitudes and the two-loop corrections to the form factor which are 
available in the massive case and it has the same structure of infrared divergences as the full 
amplitude. 
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4 Two-loop corrections to the massive Bhabha cross 
section 



Once the result for the auxiliary amplitude in Eq. (fTTjl is known, the problem is to evaluate 
the difference 



5A 
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which matches the auxiliary and the full amplitudes. We should note that though the 
different terms on the right hand side of Eq. (|32[) taken separately are infrared divergent, 
their sum can be transformed into convergent Feynman integrals. In similar way the pinch 
singularities disappear at the level of Feynman integrals in the two-loop corrections to the 
static potential after the proper infrared subtraction [35] • Thus, in principle, Eq. does 
not need to be regularized. However, it is simpler to take the available results for the 
different terms of Eq. (|3*2*|) in dimensional regularization and then to take the limit d — > 4 in 
the sum. The explicit result for the matrix elements of the two-loop and the tree amplitudes 
in dimensional regularization can be found in Ref . jH] . To disentangle the infrared divergences 
the authors of Ref. jH] used the formula suggested by S. Catani in Ref. |%B] . In the next 
section we describe how this result can be matched to Eq. (J32j) . 



4.1 Two-loop infrared matching 

In the case of pure photonic corrections the Catani formula for the structure of the infrared 
divergences takes the following form 
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which are diagonal in the chiral basis. The form of the nonsingular terms in Eq. (|3"4*|) is a 
matter of convention. We use the one of Refs. |8|4fi| 1 . It is easy to check that the above 
expression is in full agreement with Eq. (}3*Tj) . Indeed, Eq. (}3*3*j) is invariant under a redefinition 
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1 Our normalization of the operators differs from [BH5j by the overall factor 2 per loop. 
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where G and F are the two- component functions of x and e which are regular at e = 0. By 
choosing 

jP = 2/ (2)_ (/ (l) ) 2_ i? (2) (36) 



G - A (1) 

U — ^fin J 



we reproduce the structure of Eqs. (J3T1 |32J) with 

4(1) = ; 2 / (2) - (/ {1) ) 2 = if /(2) 



M< 2 > = A' (2) 



fin 



(37) 



Thus Eq. can be considered as a direct consequence of the evolution equations ([TH | I?7 ]) . 
Moreover, by analysing the evolution equations we can predict the form of the operator 
for the scattering amplitude which was not determined in Ref. I3H] but instead has been 
found by explicit calculation [S]. In fact the same is true for the four-quark amplitude where 
the Catani formula is a direct concequence of the non-Abelian evolution equations [HE] . 

By using Eq. (|3fi|) one can transform the result of Ref. [B] for the matrix element of 
Aftl into the one of 5A^. However, we prefer to use directly the result of [Hj for the 
finite part of the two-loop corrections and find the expressions for the operators I^ 1 ' and 
jj (2) corresponding to the mass regularization of the infrared divergences. To perform this 
infrared matching let us first note that the primed operators defined through Eqs. ()351 l3*Bj) 
are given by the sum of Feynman integrals corresponding to the one-loop correction to 
the amplitude and two-loop correction to the logarithm of the form factor, respectively. 
Therefore, in contrast to the original definition ()34)1 . it is straightforward to obtain the 
variation of the primed operators with the change of the infrared regularization. For 1'^ it 
consists in replacing the one-loop massless result by the massive one. For H'^ the matching 
term which relates the dimensionally regularized and the massive result is given by twice 
the difference of the nonlogarithmic terms of Eq. (|2~T|) and Eq. (|2T)|) . where the (f^) 2 /2 
contribution is subtracted. As we have already pointed out, the finite part A'^l = 8A^> 
does not depend on the regularization. Now we can perform the inverse transformation to 
the operators 1^ and H^ 2 \ Note that we are interested in the limit d = 4 and only need 
the value of the functions F and G at e = 0. The finite part of the two- loop correction to 
the amplitude A^J does not change after these transformations while 1^ and become 
logarithmic functions of electron and photon masses 
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By plugging Eq. (|38|) into Eq. (J3*3*j) we reproduce the known result for the logarithmic cor- 
rections to the amplitudes. Beside the logarithmic correction the operators (|3~%j) produce a 
nonlogarithmic contribution in Eq. ()33j) which can be considered as the matching term be- 
tween the dimensionally regularized and the massive result for the amplitudes. Note that in 
Ref. jH] the explicit result is given only for the matrix elements of the two-loop and the tree 
amplitudes rather than the expressions for the chiral amplitudes. This, however, is sufficient 
for the matching because the operators Eq. ()34|) are diagonal in the chiral basis. 



4.2 The result 

Now we are in a position to derive the result for the second order correction to the cross 
section of the massive Bhabha scattering. It can be split into three parts: 

(i) the corrections involving the soft real emission; 

(ii) the interference of the one-loop corrections to the amplitudes; 

(iii) the interference of the two-loop corrections and the tree amplitudes. 

The soft photon emission is known to factorize and the corresponding second order correc- 
tions to the cross section introduced in Eq. (jHJ) are of the following form 
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With the known one-loop corrections to the chiral amplitudes at hand (see e.g. Ref. 07j) it 
is straightforward to obtain the corresponding interference term in the cross section, Eq. (J4~5j) 
of the Appendix. The derivation of the contribution (iii) has been described in the previous 
section. Collecting all the contributions we obtain the result for the nonlogarithmic photonic 
correction which is given by Eq. ()4fij) of the Appendix. 

In the limit of small scattering angles the virtual corrections to the cross section are 
completely determined by the corrections to the electron and positron form factors in the 
t-channel amplitude [20] • We check that our result for the the virtual corrections in the the 
limit x — > reduces to 
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where ff^ are given by Eqs. (fHfl l2*Tj) with Q 2 = xs. This agrees with the asymptotic small 
angle expression given in [21] which is quite a nontrivial check of our result. The two- loop 
nonlogarithmic corrections to the cross section in the small angle limit becomes 
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Figure 1: (a) Logarithmically enhanced (dashed line) and nonlogarithmic (solid line) second order 
corrections to the differential cross section of the small angle Bhabha scattering as functions of the 
scattering angle for = 100 GeV and ln(e cut /e) = 0, in permill. (6) The same as (a) but for the 
large angle Bhabha scattering and yfs = 1 GeV. 



(a) 



where 

27 17 8 
C = — + — vr 2 - 9C(3) - ^vr 4 - 2vr 2 ln(2) = -0.744199 .... (42) 
2 8 45 

Note that our result is not valid for very small scattering angles corresponding to x,< m 2 / s and 
for almost backward scattering corresponding to 1 — x<ml/s, where the power-suppressed 
terms of the form m 2 e /t and m^/u become important. 

Our result should be combined with the Monte-Carlo evaluation of the hard brems- 
strahlung. As it has already been mentioned in many practical realizations of the Monte- 
Carlo event generators the cancellation of infrared divergences between virtual and soft real 
corrections is implemented to high orders in perturbation theory by using the exponentiation 
property discussed above. In this case the Monte-Carlo result for the cross section already 
includes a part of the second order virtual and soft real corrections. This part depends on the 
specific realization of the Monte-Carlo algorithm and should be subtracted from the result 
of the paper. 



5 Numerical estimates and summary 

Let us now discuss the phenomenological relevance of our result. The determination of the 
luminosity for the GigaZ option of ILC is most demanding with respect to the theoretical 
predictions for the small angle Bhabha scattering. It requires the accuracy at the level of 0.1 
permill [3] . At the same time the low-energy experiments aimed at the determination of the 
hadronic vacuum polarization contribution through the measurement of a(e + e~ — *> hadrons) 
require about one permill accuracy of the large angle Bhabha cross section. Such a high accu- 
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racy is necessary to reduce theoretical uncertainty due to the hadronic vacuum polarization 
contribution to the muon anomalous magnetic moment and to the value of the QED coupling 
constant at Z peak (see e.g. Ref. jlEJ)- Neither of the existing Monte-Carlo event generators 
for small |l|2|22j and large angle |5fl0|llj Bhabha scattering, which as yet do not incorporate 
the complete second order QED corrections, can guarantee the required precision. 

The second order photonic corrections to the differential cross section (a/7r) 2 dcr^ 2 ^/d(j^ 
are plotted as functions of the scattering angle for the small angle Bhabha scattering at 
y/s = 100 GeV on Fig. (JT^) and for the large angle Bhabha scattering at ^fs = 1 GeV 
on Fig. We separate the logarithmically enhanced corrections given by the first two 

terms of Eq. (j^J) and the nonlogarithmic contribution given by the last term of this equa- 
tion. All the terms involving a power of the logarithm \n(e cu t/e) are excluded from the 
numerical estimates because the corresponding contribution critically depends on the event 
selection algorithm and cannot be unambiguously estimated without imposing specific cuts 
on the photon bremsstrahlung. The actual impact of the two-loop virtual corrections on 
the theoretical predictions can be determined only after the result of the paper is consis- 
tently implemented into the Monte-Carlo event generators. Nevertheless, the above naive 
procedure can be used to get a rough estimate of the magnitude and the structure of the 
corrections. We observe that for scattering angles 6<,18° and 6^166° the nonlogarithmic 
contribution exceeds a benchmark of 0.5 permill which makes it relevant for the luminosity 
determination at the low-energy electron-positron colliders. For the small scattering angles 
the second order correction reaches a few permill in magnitude. Here we should note that 
BHLUMI event generator [22] used for luminosity determination at LEP includes the total 
leading logarithmic second order contribution enhanced by the factor \r?(t/m 2 e ) as well as the 
bulk of the subleading contribution. In fact the remaining part of the subleading photonic 
corrections to the cross section has been computed in the small-angle approximation [0] but 
has not been included in the code. According to Ref. || this missing correction amounts for 
approximately 0.14 permill for the energy and scattering angles characteristic to LEP, which 
is relevant for the GigaZ accuracy. 

To get the total second order correction (without the hard bremsstrahlung) our result 
should be combined with the fermionic contribution. In Fig. (j2J) we plot the second order 
photonic contribution against the fermionic one in the case of one light flavor. The fermionic 
contribution incorporates the second order corrections with one closed fermion loop including 
the single soft photon emission ^3] and the contribution due to the emission of the soft real 
electron-positron pair of the energy below a cutoff e e ^ -C s. The latter has been computed 
in Ref. [49J in the logarithmic approximation and cancels the artificial In 3 [m 2 e / s) term of the 
closed fermion loop contribution. Note that we do not include a trivial contribution with two 
closed fermion loops which can be taken into account through the one-loop renormalization 
group running of a in the tree amplitudes. 

Finally, we would like to mention the electroweak corrections to Bhabha scattering which 
can be important at the considered level of accuracy. The one-loop correction is well known 
|47j . For the large angle scattering above the electroweak scale, however, the two- loop 
electroweak corrections could be important. In the case of e + e~ — * fi + ^~ annihilation the 
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Figure 2: (a) Photonic (solid line) and fermionic (dashed line) second order corrections to the 
differential cross section of the small angle Bhabha scattering as functions of the scattering angle 
for yfs = 100 GeV and ln(e cu t/e) = ln(e^ e je) =0, in permill. (b) The same as (a) but for the 
large angle Bhabha scattering and y/s = 1 GeV. 



corrections enhanced in the high energy limit by a power of the large logarithm ln(M 2 /s), 
where M stands for W or Z boson mass, have been computed in |37|5U|51|52j . They dominate 
the electroweak corrections for the energies a/s>,500 GeV characteristic to ILC. Due to the 
strong numerical cancellations between the terms with different powers of the large logarithm 
the total electroweak logarithmic two-loop contribution does not exceed a few permill in this 
energy region. This analysis can be generalized to the large angle Bhabha scattering by 
adding the t-channel contribution. 

To conclude, we have derived the two-loop radiative photonic corrections to Bhabha 
scattering in the leading order of the small electron mass expansion up to nonlogarithmic 
term. Together with the result of Ref. for the fermion loop corrections our result 

gives a complete expression for the two-loop virtual corrections. It should be incorporated 
into the Monte Carlo event generators to match the demands of the present and future 
electron-positron colliders for the accuracy of the luminosity determination. 
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Appendix 

The one-loop virtual photonic correction to the normalized cross section reads 
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(43) 



where f(x) is given by Eq. (J7J). The correction due to the single soft photon emission reads 
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(44) 



The second order contribution to the cross section due to the interference of the one-loop 
virtual corrections to the amplitude reads 
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where the trivial terms proportional to ln(A 2 /m 2 ) and ln(m 2 /s) are omitted. 
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The total second order nonlogarithmic photonic contribution to the normalized cross 
section reads 
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where £ e = [1 - In (x/(l - x))) In (e cut /e). 
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